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Abstract 

The calculation of the electronic conductance of nanoscale junctions from first principles is a long-standing problem in the field of 
charge transport. Here we demonstrate excellent agreement with experiments for the transport properties of the gold/alkanediamine 
benchmark system when electron-electron interactions are described by the many-body GW approximation. The conductance 
follows an exponential length dependence: G„ = exp(-pw). The main difference from standard density functional theory (DFT) 
calculations is a significant reduction of the contact conductance, G^, due to an improved alignment of the molecular energy levels 
with the metal Fermi energy. The molecular orbitals involved in the tunneling process comprise states delocalized over the carbon 
backbone and states localized on the amine end groups. We find that dynamic screening effects renormalize the two types of states 
in qualitatively different ways when the molecule is inserted in the junction. Consequently, the GW transport results cannot be 
mimicked by DFT calculations employing a simple scissors operator. 



Introduction 

The conductance of a molecule sandwiched between metallic 
electrodes is sensitive to the chemical and electronic structure of 
the molecule as well as the detailed atomic structure of the 
metal-molecule contact. Variations in the contact geometry 
beyond experimental control lead to an undesired spread in the 



measured conductance properties. For the most commonly used 
anchoring group, -thiol, these effects are rather pronounced due 
to the many possible contact geometries resulting from the 
strong Au-S interaction. Amine groups have been shown to 
produce more well-defined transport properties [1], which can 
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be understood from the relatively weak AU-NH2 bond leading 
to larger structural selectivity [2]. 

Even for a given junction geometry, a quantitatively accurate 
description of electron transport from first principles remains a 
formidable task. Numerous studies based on density functional 
theory (DFT) have shown a significant overestimation of 
conductance relative to experimental values [3-14] (an 
exception to this trend occurs for small molecules, such as H2 
[15-17] and CO [18,19], which chemisorb strongly to the 
electrodes, resulting in resonant transport through broad, 
partially filled resonances). The inability of DFT to 
describe off-resonant tunneling in the simplest molecular 
junctions limits the predictive power of the DFT-based 
approach to qualitative trends. It is now broadly accepted that 
the failure of DFT is mainly due to its incorrect description of 
the molecular energy levels. Indeed, physically motivated 
correction schemes have shown that much-improved agreement 
with experiments can be obtained after shifting the DFT 
molecular energy levels [13,14]. Such corrections are supposed 
to remove the self-interaction errors inherent in standard DFT 
exchange-correlation (xc) functionals [20-22] and account for 
image charge effects induced by the metal contacts. The draw- 
back of the approach is that it assumes a weak coupling 
between molecular orbitals and metal states and treats the 
image-plane position as a free parameter. 

The (self-consistent) GW approximation [23], which is rooted 
in many-body perturbation theory, was recently found to yield a 
considerable improvement over DFT for the conductance of 
gold/benzenediamine junctions [24]. Physically, the GW 
approximation corresponds to Hartree-Fock theory with the 
bare Coulomb interaction v = l/|r - r'| replaced by a dynami- 
cally screened Coulomb interaction ^(co) = 8~^(co)v. In contrast 
to standard DFT, the GW approximation is almost self-inter- 
action free [25] and includes screening effects through the 
correlation part of the self-energy [26-28]. As a consequence, it 
provides quantitatively accurate predictions of energy gaps in 
systems with highly diverse screening properties, ranging from 
isolated molecules [29,30], through to semiconductors [31] and 
metals [32]. The broad applicability of the GW approach 
becomes particularly important for a metal-molecule interface 
where the electronic structure changes from insulating to 
metallic over a few angstroms. 

In this work we use the GW approximation to study the role of 
exchange-correlation effects for the energy-level alignment and 
electron transport in short alkane chains coupled to gold 
electrodes through amine linker groups. The gold/alkane 
junction is a benchmark system for molecular charge transport 
and has been exhaustively investigated experimentally [1,12,33- 



44]. We focus here on the amine-linked alkanes to avoid the 
uncertainties related to the gold-thiol contact geometry, which 
is presently under debate [45-50]. We note that very recently it 
was shown that alkanes can be bound directly to gold electrodes 
without the use of anchoring groups [51]. 

The transport mechanism in (short) saturated molecular wires is 
coherent tunneling through molecular orbitals with energy far 
from the Fermi energy. The trend of conductance versus chain 
length (n) thus follows an exponential law of the form 

G„ = exp(-p«) (1) 

Recent experimentally reported values for the decay constant P 
of alkane-a,ca-diamine/gold junctions are in the range 0.9-1.0 
per C atom [1,12,36], but earlier measurements also showed 
values around 0.8 [35]. Although previous studies based on 
DFT yielded P values within the experimental range, the contact 
conductance, G^, is typically overestimated by around an order 
of magnitude [3-10]. A study based on the many-body configur- 
ation interaction method has shown similar p values, but 
slightly reduced values, as compared to DFT [52]. By 
comparing DFT and GW calculations for C„-alkanediamine 
molecules with « = 2,4,6 we show that the erroneous Gc values 
are a result of the incorrect level alignment in the DFT calcula- 
tions. Indeed, GW yields a G^ in close agreement with the 
experimental values. We find a pronounced orbital and length 
dependence of the quasiparticle (QP) corrections to the DFT 
energies, resulting from the different shape and localization of 
the molecular orbitals. The QP corrections range from -0.5 to 
-2.5 eV and can be qualitatively explained from a classical 
image-charge model. 

Method 

The junction geometries were optimized by means of the real- 
space projector- augmented wave method GPAW [53,54] with a 
grid spacing of 0.2 A and the PBE functional for exchange and 
correlation (xc) [55]. The molecules were attached to Au(l 11) 
surfaces, modeled by an eight-layer-thick 4x4 slab, through 
small four-atom tips as shown in Figure la. The surface Bril- 
louin zone was sampled on a 4 x 4 Monkhorst pack ^-point 
grid, and the structures including molecule, Au tips, and outer- 
most Au surface layers were relaxed until the residual force was 
below 0.03 eV/A. We considered w-alkanediamine junctions 
with /7 = 2, 4 and 6. The key structural parameters can be found 
here [56]. For calculations of the molecules in the gas-phase, we 
include 16 A of vacuum between molecules in the repeated 
supercells. All transport calculations where performed 
according to the method described in detail in [24]. In brief, we 
employ a basis set of numerical atomic orbitals corresponding 
to double-zeta plus polarization (DZP) for the Au atoms and 
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double-zeta (DZ) for the atoms of the molecules. We use rather 
diffuse basis functions with a confinement-energy shift of 
0.01 eV. This ensures that the calculated work function of 
Au(lll) and the Kohn-Sham energy levels of the molecular 
junction are within 0.1 eV of those obtained from accurate grid 
calculations [24]. The transmission function is obtained from 
the Meir-Wingreen transmission formula [57,58] 

Tie) = Tv[G' (s)r^ (£)G'' (£)r^ (s)] . (2) 

The retarded Green's function of the extended molecule is 
calculated from 

G' (E) = [(E + iri)S -Hq+ V,, -AV^ [G] 
-Si(£)-2«(£)-2,JG](£)rl 

Here S, Hq, and V^c are the overlap matrix, Kohn-Sham Hamil- 
tonian and the local xc-potential in the atomic orbital basis, res- 
pectively; Tj is a positive infinitesimal. 




Figure 1: (a) Supercell used to model the gold/alkanediamine 
junctions. Similar supercells were used for n = 2 and n = 4 (not shown) 
and the key structural parameters can be found here [50]. The black 
box indicates the region of the extended molecule within which the GW 
self-energy is evaluated, (b) Isosurface of the electron-density differ- 
ence between a DFT ground-state calculation and a constrained DFT 
calculation with one less electron on the molecule. 

The lead self-energy, ^l/r^ incorporates the coupling to the left 
and right electrodes and is obtained by standard techniques [59]. 
The term AF// gives the change in the Hartree potential relative 
to the DFT Hartree potential already contained in Hq. Finally, 
the last term is the many-body xc self-energy, which in this 
work is either the bare exchange potential, V^, corresponding to 
Hartree-Fock, or the GW self-energy. As indicated both the 
Hartree potential and the xc self-energy depend on the Green's 
function. The latter is evaluated fully self-consistently using a 
simple linear mixing of the Green functions. We represent all 



energy dependent quantities in Equation 3 on a large energy 
grid ranging from -200 eV to 200 eV with an energy-grid 
spacing of 0.01 eV. 

The GW self-energy is evaluated for the extended molecule 
(indicated by the box in Figure la). However, only the part 
corresponding to the molecule is used while the remaining part 
is replaced by the DFT xc-potential. This is done to include 
nonlocal correlation (image-charge) effects from the electrodes 
in the GW self-energy of the molecule while preserving a 
consistent description of all metal atoms at the DFT level. We 
have verified that the calculations are converged with respect to 
the size of the extended molecules, see [24] for more details. 
The basis functions on the Au tip atoms extend a fair distance 
into the Au electrode and thus also describe screening effects to 
some extent in this region. The rather fast convergence of the 
screening effects with extended molecule size can be probed 
directly through a simple approach based on a constrained DFT 
calculation where the number of electrons on the molecule is 
constrained to be one less, or one more, than the number of 
electrons in a ground state DFT calculation. We measure the 
number of electrons on the molecule using a Mulliken charge 
analysis. The constrained DFT calculations are performed as 
self-consistent DFT calculations with the rigid shift of the 
molecular orbitals adjusted until the number of electrons is one 
less, thereby probing the highest occupied molecular orbital. 
We show in Figure lb the isosurface of the electron-density 
difference obtained from a ground state DFT calculation and a 
constrained DFT calculation. The change in the electron density 
in the metal can be seen as the formation of an image charge 
(red color), which to a large extent is localized on the Au tip 
atoms closest to the molecule. 

Results and Discussion 
Energy-level alignment 

The alignment of the molecular energy levels relative to the 
electrode Fermi level is of great importance for the transport 
properties of molecular junctions and seems to be the dominant 
effect at low bias voltage. At higher bias voltages, many-body 
calculations on small model systems suggest that electron corre- 
lations induce additional shifting and broadening of the 
molecular levels, which can also affect the transport properties 
[60]. Here we focus on the low-bias regime and postpone 
consideration of the finite-bias effects to a later study. 

The molecular orbitals (MOs) of the alkanediamine chains 
comprise states that are delocalized over the carbon backbone 
and states that are localized on the NH2 end group. We shall 
consider the highest occupied molecular orbital (HOMO) and 
HOMO-2 as representatives for the two classes of states 
(Figure 2a). We note that the HOMO-1 is similar to the 
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HOMO, with slightly lower energy given by the coupling of the 
two end groups across the wire. In Table 1 and Table 2 we list 
the energy of the HOMO and HOMO-2 calculated with DFT- 
PBE, Hartree-Fock (HF), and GW for the molecules in the gas- 
phase and in the junction. 
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Figure 2: (a) Isosurfaces for the HOMO and HOMO-2 orbitals of the 
Ce-alkanediamine molecule, (b) HOMO and HOMO-2 PDOS in the 
junction obtained from DFT-PBE (top), GW (middle) and DFT with a 
generalized scissors operator (bottom), (c) Quasiparticle corrections to 
the Hartree-Fock levels in the gas-phase (open symbols) and in the 
junction (filled symbols) as a function of molecular length n. The 
HOMO and HOMO-2 are denoted by circles and squares, respective- 
ly, (d) Same as (c) but for DFT rather than Hartree-Fock. 



In the gas-phase, all three methods predict the HOMO energy to 
be almost independent of molecular length. This is clearly due 
to its end-group-localized character. In contrast, the energy of 
the HOMO-2 level shifts upward in energy as the molecular 



Table 1: Calculated HOMO and HOMO-2 energies aligned to the 
vacuum level and in units of eV. 



method 



DFT-PBE 



HF 



GW 



orbital 




A7 = 2 


n = 4 


n = 6 


HOMO 




-4.9 


-5.1 


-5.1 


HOMO- 


-2 


-8.5 


-8.2 


-8.0 


HOMO 




-10.2 


-10.5 


-10.5 


HOMO- 


-2 


-13.3 


-12.9 


-12.8 


HOMO 




-8.5 


-8.6 


-8.6 


HOMO- 


-2 


-12.2 


-11.8 


-11.6 



Table 2: Calculated HOMO and HOMO-2 energies in the junction 
relative to the electrode Fermi level. 



method 


orbital 


A? = 2 


A? = 4 


n = 6 


DFT-PBE 


HOMO 


-4.3 


-4.2 


-4.4 




HOMO-2 


-6.5 


-5.9 


-5.7 


HF 


HOMO 


-8.1 


-8.1 


-8.3 




HOMO-2 


-10.4 


-9.9 


-9.7 


GW 


HOMO 


-4.8 


-4.9 


-5.2 




HOMO-2 


-8.4 


-8.2 


-8.2 



length increases. This reflects its extended nature and can be 
interpreted as a band-discretization effect. To the best of our 
knowledge no experimental results exist for the ionization 
potential of alkanediamine molecules. However, for the closely 
related butane molecule (C2 alkane with CH3 end groups) we 
obtain a GW-calculated HOMO energy of -11.4 eV in very 
good agreement with the experimental ionization potential of 
11.2 eV [61]. In comparison, the DFT-PBE HOMO energy is 
severely underestimated at -7.9 eV. This finding agrees well 
with previous studies on a broader range of small molecules 
[24,29,30]. 

In the junction, the molecular orbitals, | were obtained by 
diagonalizing the DFT Hamiltonian corresponding to the mole- 
cule. The projected density of states (PDOS) of such a state is 
then given by the spectral function, -1 / 7T;Im(\|/^ | G^(E) | 
where G is the appropriate Green's function (calculated with 
DFT, HF, or GW). The level position is defined as the first 
moment of the PDOS. Figure 2b shows the PDOS for the 
HOMO and HOMO-2 for the Cg-alkanediamine junction as 
calculated with DFT-PBE (upper panel) and GW (middle 
panel). The lower panel shows the PDOS obtained from a DFT 
calculation where the molecular levels have been shifted to 
match the GW levels, i.e., after adding to the Kohn-Sham 
Hamiltonian a generalized scissors operator of the form 
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n 

Here, the s^^ denote the QP energy obtained from the GW 
calculation. We see that the main features of the GW spectral 
function can be well reproduced by the shifted DFT Hamil- 
tonian, although small differences remain. A similar conclusion 
was reached in [24] for a gold/benzenediamine junction. 

The molecular orbital energies from a GW calculation include 
the dynamical response of the electron system to an added or 
removed electron through the correlation part of the self-energy. 
In general, correlations tend to shift the filled levels upwards 
and the empty levels downwards relative to the bare 
Hartree-Fock energies. This is because the inclusion of 
screening reduces the energy cost of removing/adding electrons 
to the molecule. When a molecule is brought into contact with a 
metallic junction its environment changes from insulating to 
metallic. This implies extra screening of an added or removed 
electron, which will cause the filled levels to shift upwards and 
the empty levels to shift downwards even more than for the 
isolated molecules, i.e., the gap will shrink relative to its gas- 
phase value. It has been shown previously that DFT in 
(semi)local approximations and Hartree-Fock completely miss 
this important effect [26-28]. 

In Figure 2c and Figure 2d, we show the QP corrections to the 
HF and DFT Kohn-Sham energy levels as a function of 
molecular length. The results for the HOMO and the HOMO-2 
are denoted by circles and squares. We notice first that the QP 
corrections are very significant with absolute values reaching 
almost 4 eV and with a pronounced orbital and length depend- 
ence. The Hartree-Fock QP corrections are all positive showing 
that HF places the occupied levels lower than predicted by GW. 
This is in contrast to the corrections to the DFT levels, which 
are all negative, in agreement with the well-known underesti- 
mation of ionization potentials as predicted from the negative 
Kohn-Sham HOMO energy obtained using LDA or GGA func- 
tionals. In contrast to Hartree-Fock the Kohn-Sham QP correc- 
tions are smaller for molecules in the junction compared to the 
gas-phase. In fact, the position of the HOMO level by DFT is 
relatively close the GW level position and only lies 0.5-0.8 eV 
higher. The fact that the DFT-PBE description of molecular 
energy levels is much better in the junction than in the gas- 
phase agrees with previous findings [24,27,62] and can be 
explained from the origin of the PBE functional in the homoge- 
neous electron gas [63]. 

It is instructive to consider the shift in the molecular energy 
levels due to correlation effects coming from the metal 
electrodes. In simple terms this corresponds to the shift induced 



by image-charge effects. In order to isolate the part of the corre- 
lation energy originating from the metallic electrodes we define 
the quantity 

^^corr,^ V*^^ /junction V^^ /gas-pnase v^.' 



which is shown in Figure 3. 
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Figure 3: Change in the correlation energy of the HOMO and 
HOMO-2 energy levels when the molecules are taken from the gas- 
phase into the junction. This value represents the shift of molecular 
levels due to the enhanced screening provided by the metallic 
electrodes (an image-charge effect in simple terms). 



The result can be understood qualitatively by considering a 
classical model to account for the screening effect of the 
electrodes. Classically a charge distribution close to any surface 
will experience an image potential. The strength of the image 
potential in general depends on the dielectric constant of the 
surface material and the local geometrical shape of the surface. 
Here we model the Au electrodes as perfect metals. The image 
potential for a point charge halfway between two metal surfaces 
separated by a distance L is ~ 10.0/L (eV-A) [64]. This predicts 
that the image-charge effect is proportional to \IL. The Au tip 
atoms in our simulations are about 8 A apart for the C2-alkane- 
diamine junction giving a rough estimate of 1.3 eV for the 
image-charge effect in qualitative agreement with the GW 
calculations. The HOMO experiences a larger image-charge 
effect than the HOMO-2, which can be understood from the 
fact that its charge density is located closer to the metallic 
surfaces. In the limit of an infinitely long wire the HOMO-2 
will be spread out over the entire molecule and the image- 
charge effect should vanish. On the other hand, in this limit the 
HOMO would stay localized near the surface and therefore 
approach a nonzero constant image charge potential. If we 
model the HOMO charge density as a point charge of half an 
electron on each of the amine groups we can estimate this 
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limiting constant to be 3.6/(2d) (eV-A), where d is the distance 
to the nearest metal surface. Taking d to be about the same 
length as the Au-N bond (2.34 A) gives a limiting value esti- 
mate of 0.8 eV. Again, this seems to be in qualitative agree- 
ment with our GW findings. 

Finally, we discuss the coverage dependence of the energy level 
position for alkanediamine-Au junctions. It was shown in [65] 
that the DFT level position for amine-linked molecules is 
strongly dependent on the coverage. In contrast to the screening 
(image-charge) effects as discussed above, which appear in the 
correlation part of the self-energy, this is a purely electrostatic 
effect resulting from the localized surface dipoles formed at the 
AU-NH2 bond. To investigate the dependence of the energy 
levels on the coverage for our junctions, we have performed 
DFT calculations for a range of transverse supercell dimen- 
sions for the geometry shown in Figure 1 . In Figure 4 we show 
the PDOS of a methylene unit in the central part of the mole- 
cule for transverse supercells with 2 x 2 up to 8 x 8 surface 
atoms. 

The PDOS peaks and band edge shift up in energy as the trans- 
verse supercell size is increased, in agreement with the results 
of [65]. The inset shows the energy shift obtained from reading 
off the shift in the PDOS as a function of coverage defined as 
Tj = l/A/'surface? where A/'surface is the number of surface atoms. 
For a supercell of size 3x3 and larger, the shift is seen to be 
directly proportional to the coverage, as expected for a two- 
dimensional array of dipoles [66]. This allow us to extract the 
electrostatic shift corresponding to the single-molecule limit. 
We find that the electrostatic energy shift when going from a 
4x4 supercell to the single-molecule limit is indeed significant 
with a value around 1 eV. 



Transport calculations 

The transmission functions of the C2-, C4- and Cg-alkanedi- 
amine junction geometries were calculated by using the GW 
and the PBE xc potential as approximations for in 
Equation 3. To include the coverage dependence, we simulated 
the low-coverage limit rj = 0^ by performing calculations for the 
4x4 junction (corresponding to rj = 1/16) with all molecular 
levels shifted up by 1 eV through a simple scissors-operator 
self-energy. 

The transmission function calculated by GW for a coverage of 
Tj = 1/16 is shown in Figure 5 on a logarithmic scale. The trans- 
mission functions for different molecular lengths have very 
similar shapes in the important region near the Fermi level E^, 
however, the magnitude is increasingly suppressed as a func- 
tion of the molecular length. The similarity of the transmission 
functions may at first seem surprising since we have shown that 
the position of the molecular energy levels shows some length 
dependence. In particular the HOMO level was found to 
decrease in energy by 0.5 eV when n increases from 2 to 6 
(Table 2). This shift is indeed visible in the transmission func- 
tion in the range -4.0 to -6.0 eV where the HOMO is located. 
On the other hand the features in the transmission function 
around the Fermi level are determined by the local electronic 
structure of the Au tips. 

The zero-bias conductance is obtained from the transmission 
function at the Fermi level, G = G()T(Ey) where Gq = 2e^/h is 
the unit of quantum conductance. The zero-bias conductance is 
plotted in Figure 6 as a function of the molecular length. We 
have also included the DFT results for comparison. The dashed 
lines show the best fits to the exponential form G„ = 
Gcexp(-p«). 
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Figure 4: The molecular valence-band edge (or HOMO) as a function 
of coverage for an n = 6 alkane. The numbers indicate the number of 
surface atoms. 




Figure 5: The transmission function calculated by GW for a molecular 
length of n = 2, n = 4 and n = 6 at a coverage of n = 1/16. 
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Figure 6: Calculated conductance plotted as a function of the 
molecular length for a coverage corresponding to 4 x 4 Au atoms per 
molecule (n = 1/16) and extrapolated to the single-molecule limit 
(n = 0). 



The values for p and corresponding to the single-molecule 
limit (ri = 0) and 4x4 Au atoms per molecule (r\ = 1/16) are 
reported in Table 3 together with the experimental values. We 
note that the contact conductance was not stated in [12], but has 
been estimated by extrapolation to zero molecular length from 
the reported experimental data. When comparing to experiment 
it should be kept in mind that experiments are often performed 
in solution at room temperature and are subject to variations in 
the detailed atomic structure. However, it has been shown that 
amine-linked molecules bind preferentially to under-coordin- 
ated Au atoms, such as in the structures considered here, and 
show a relatively narrow conductance distribution [1]. 



Table 3: Calculated contact conductance Gc in units of Gq and expo- 
nential decay constant (3 per carbon atom; n denotes the coverage. 




DFT 


GW 


exp. 


n 


1/16 0 


1/16 0 




Go 


0.11 0.10 


0.02 0.04 


0.02^, 
0.030^ 


p 


0.99 0.92 


1.01 1.02 


0.97^, 
0.93^ 
0.91^ 


"[12] (Gc 
^[1]. 


is not stated, but is estimated from the reported data); ^ [31]; 



The rather weak effect of coverage on the conduction prop- 
erties is in agreement with the findings reported in [67] where a 
C4-alkanediamine/Au junction with 3x3 and 4x4 surfaces 
was considered. Our DFT results are in reasonable agreement 
with previous DFT studies showing decay factors in the range 
0.83-1.01 and conductance resistances in the range (0.09-0.28) 
Go [8-10]. While the p values obtained with GW are rather 



close to the DFT calculated ones, the contact conductance is 
reduced by a factor of 3-5, depending on the coverage. This is a 
direct result of the molecular levels lying further away from 
(by 0.5-2.5 eV, Figure 2d) in GW compared to DFT. 

Conclusion 

We have unraveled the important role of exchange-correlation 
effects for the energy-level alignment and low-bias conduc- 
tance of gold/alkanediamine molecular junctions. Based on 
many-body GW calculations we found that the origin of the 
overestimation of the contact conductance, G^, by standard DFT 
is due to the incorrect energy-level alignment in the junction. 
The absence of self-interaction and the inclusion of image- 
charge screening effect through the GW self-energy improves 
the description of the energy levels and yields values for G^ and 
the decay constant p that are in good agreement with experi- 
mental results. The quasiparticle corrections to the DFT energy 
levels showed a significant orbital dependence ranging from 
-0.5 eV to -2.5 eV due to the different shape and localization 
of the molecular orbitals. Our results demonstrate that quantita- 
tively accurate calculations of conductance from first-principles 
are feasible, although computationally demanding. 
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